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We consider the densities of clusters, at the percolation point of a two-dimensional system, which 
are anchored in various ways to an edge. These quantities are calculated by use of conformal field 
theory and computer simulations. We find that they are given by simple functions of the potentials 
of 2-D electrostatic dipoles, and that a kind of superposition cum factorization applies. Our results 
broaden this connection, already known from previous studies, and we present evidence that it is 
more generally valid. An exact result similar to the Kirkwood superposition approximation emerges. 

PACS numbers: 64.60.Ak,68.35.Rh,68.18. Jk 



Percolation in two-dimensional systems is an area with 
a long history, which remains under very active current 
study. A plethora of methods has been applied to critical 
2-D percolation, including conformal field theory (CFT) 
0, modular forms 0, computer simulation other 
field-theoretic methods Stochastic Lowner Evolution 
(SLE) processes and other rigorous methods [|| . (The 
literature is so extensive that we have cited only a very 
few representative works.) 

It has long been known that the behavior of perco- 
lating systems at or near the percolation point is inti- 
mately connected with two-dimensional electrostatics, or 
more precisely the Coulomb gas (Gaussian model) of field 
theory. In particular, Coulomb gas methods have given 
many important results, for instance crossing formulas 
on a torus • Recently, Smirnov has proved that the 
crossing probability in critical percolation in an equilat- 
eral triangle may be described as the boundary value of a 
certain very simple problem in 2-D electrostatics. In this 
Letter, we deepen and extend this connection by demon- 
strating that certain correlation functions, which specify 
the density of clusters at the percolation point, are given 
as simple functions of the potentials of two-dimensional 
dipoles. 

In addition we show that a kind of superposition ap- 
plies. An exact result resembling the Kirkwood superpo- 
sition approximation emerges. 

We find indications that our results also apply when 
the anchor points are not on the edge, and hold for the 
critical Potts models. 

Our results are found via CFT [|| with verification by 
computer simulations. 

Our work considers the (average) density of anchored 
critical clusters in 2-D percolation. The clusters are in 
the upper half-plane, and are constrained to touch the 
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real axis at a specified point or points, or in an inter- 
val between two given points. By conformal invariance, 
the corresponding results for any compact region may be 
obtained by using the appropriate transformation. The 
related problem of the density of clusters in statistical 
systems in other geometries with fixed boundaries has 
been investigated in 

The first case we consider anchors the clusters at a 
single point, which we may take as the origin. To find 
their density, we use Cardy's CFT analysis Q of cross- 
ing probabilities in critical 2-D percolation for systems 
with an edge. This approach makes use of the Q-state 
Potts model and the boundary operator 4>\,2{x), which 
changes the boundary conditions from fixed to free at a 
point x on the real axis. For percolation (Q — > 1) it has 
conformal dimension hi t a = 0. Two such boundary oper- 
ators </>i 2(s ; i)</ , i, 2(^2) may be used to count the number 
of clusters touching the real axis for x\ < x < x-i [Til ]. 

If we bring X\ and X2 together, the operator product 
expansion (OPE) is applicable. Two terms arise. One 
involves the unit operator, which has conformal dimen- 
sion /ii i = 0, and the other <j)(x) = 4>i^(x). This latter 
operator has conformal dimension h$ — ft.1,3 = 1/3 > 
0, and thus (since the corresponding term vanishes in 
this limit) represents the clusters anchored between the 
points. (Since the boundary spins are fixed, in this limit 
there must be exactly one anchored cluster, as was also 
argued in The dimension 1/3 also governs the van- 

ishing of the crossing probability for, e.g., a long thin 
rectangle. The other term, from the unit operator, rep- 
resents clusters that are anchored on the real axis every- 
where except at X\ < x < x 2 . In this case the density 
becomes independent of it should. 

The operator ip(z), which measures the cluster density 
at a point z, has dimension — 5/96. The correspond- 
ing conformal operator is identified as ip = 03/2.3/2- This 
is the Q — > 1 (m — * 2) limit of the magnetization operator 
0(m±i)/2,(m+i)/2 f° r the Q-state Potts model [l3Hl4j . 

The cluster density with one anchor point is therefore 
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given by 
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where 



h/2 



and 



denote half-plane and full plane 



expectation values respectively, the second line follows 
because half-plane correlation functions are given by full- 
plane correlators with "image" operators [l4j, and the 
third line makes use of the standard form for the three- 
point function. (Here and below we ignore certain mul- 
tiplicative constants.) 

Simplifying and introducing polar coordinates gives 



Ppt(z) 
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It follows that p pt — > when z approaches the real axis 
or oo. It vanishes as p p t ~ y~ 7 / 16 when y = 3(z) — * oo, 
and as y 11 / 48 when y — > (for x ^ 0), so that the density 
has a maximum as a function of y. An exception occurs 
for the limit z — ► 0, where p p i diverges as r~ 7 / 16 . 

An infinite value for the density in a real system is 
of course not possible. Its value on any lattice, for in- 
stance, will be bounded via a cut-off involving the lattice 
constant. 

On the other hand, (JTJ may be expressed as 



P P t(z) 
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where &dip(z; xq) = 2y/ ((x — xq) 2 + y 2 ) may be taken as 
either the real (physical) part of the complex potential 
of an 2-D unit point dipole in the y direction located at 
position xq on the real axis or the imaginary part when 
the dipole is in the —x direction. 

Note the contours <&dip(z; xq) = a are circles of radius 
1/a, centered at (xo,yo = 1/a) and all tangent to the 
real axis at xq. Hence the contours of constant y 5 ^ 48 p p t 
are likewise circles. The contours of p p t itself have the 
form of circles deformed by pushing against the real axis. 
Thus lines parallel to the y-axis pass "over the shoulder" 
of the density and have a maximum (when x ^ 0). 

The above results are a nice example of the power of 
conformal invariance. Up to a multiplicative constant, 
the density p p t(z) is completely determined by only two 
numbers, the conformal dimensions 1/3 and 5/96. 

Functions of the same form as Q arise generally in 
the related problem of the evaluation of two-point cor- 
relation functions of a single operator in the upper half- 
plane when one of the operators approaches the real axis 

Comparison of J3J) with computer simulations is excel- 
lent, as described in detail below. 



It is interesting to note that 10, when transformed via 
a conformal map w — w(z), preserves its form. It remains 
a product of two factors. The first is a function of w only, 
the form of which depends on the size and shape of the 
new region, but is independent of the anchoring point 
w(0). This is multiplied by &dip(w; u^O)) 1 / 3 , the (real or 
imaginary) potential of a dipole at the mapped anchor 
point w(0) in the new region, with grounded boundary. 

The next case of interest has the cluster anchored at 
the two points x a and x b = x a + D. By the arguments 
above, the density with two anchor points is 



P2 P ts(z;x a ,x b ) = (4>(x a )4>(x b )ip{z)tp(z)) . 



(4) 



Following the standard treatment (see j^] or [HI), we 
rewrite the correlation function in (0J so that 

P2 P ts(z;x a ,x b ) = (x b -x a )- 2h *{z- z)- 2h ^F{ V ) 7 (5) 

where the cross-ratio n = [{x b — x a )(z — z)]/[(z — x a )(z — 

Xb)]. 

Since <f> = ^.3 is a level-three operator, by standard 
CFT analysis, F{rj) satisfies the third-order differential 
equation 



lOSr? 2 ^ - l) 3 -^ + 72(77 - l) 2 T7(5r? - 1)-^ 
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3(77 - l)(35ry 2 - 24)— + 577(77 - 2) 
a?7 



F(n) = 0. (6) 



Since percolation is a logarithmic CFT, the use of a 
differential equation like © , which arises from a Virasoro 
null vector, can be problematical. However this appears 
not to be the case here (see and [lUp. 

We can find the solution to @ of interest by making 
use of the OPE of (j)(x a )(j){x b ) as x a — > x b . This has 
three terms, involving either the unit operator, <f>(x a ) or 
4>i,5( x a)- On the other hand, P2pts(z',x a ,x b ) should be 
proportional to p p t{z;x a ) in this limit. The first OPE 
term clearly gives a contribution which is independent 
of x a and thus may be discarded. The operator <pi^ 
has been arg ued to create two clusters separated by a 
dual line |l2j . and thus cannot contribute to p pt . What 
remains may be expressed in the interesting factorizcd 
form 



P2 P ts(z; x a ,x b ) 



y 



5/^pT/3(®dip(z-,Xa)<f>dip(z;Xb)) 
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Clearly (7J reduces to ©, as expected, if we take D — > 
(and set x a — 0). Note also that in the limit z — ► x a or 
z — > x b , p2 P ts — 1/-D 2 / 3 , which is just the (unnormalized) 
probability of a cluster connecting x a and x b . (The latter 
limit also implies that ip reduces to <f> on a free boundary, 
as discussed further below.) 

Multiplying by y 5 / 48 and comparing J3J| and J7J| we see 
that the logarithm of the latter is linear in the logarithm 
of the former. Thus a kind of superposition applies. 
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To test these predictions, we carried out simulations on 
square lattices of sizes 255 x 255, 511 x 511, 1023 x 1023, 
and 2047 x 2047 sites using bond percolation at the crit- 
ical point p = 1/2, and keeping track of the wetted sites 
connected to the two anchor points, positioned 3/8 and 
5/8 the way up the j/-axis (e. g., in the 255 x 255 case, 
the two points separated by ±32 lattice spacings from 
the center). All boundaries were free. We used a cluster 
growth algorithm starting from either anchor point, and 
if the cluster touched the other anchor point, we also 
averaged that density separately. Thus, we found the 
density of clusters touching either one or both anchor 
points. In Fig. ^ the upper left panel shows the density 
of points touching the lower anchor point, and the upper 
right shows the density of clusters touching the upper an- 
chor point. These densities are normalized to be one at 
the anchor point. The lower left panel shows the density 
of clusters touching both anchor points, normalized simi- 
larly so that it is one at either anchor point. In the lower 
right panel, we show the square root of the product of the 
two simulated one-point densities divided by the proba- 
bility that the two anchors are connected together, finally 
multiplied by a constant C to make the contours agree 
quantitatively with those of the lower left panel. The ef- 
fective value of C is one when z is at either anchor point, 
by definition. However, it changes to C = 1.030 ± 0.001 
within a few lattice spacings, independent of lattice size. 
We interpret C below. 

We find quantitatively similar results for site percola- 
tion on square and triangular lattices, and bond percola- 
tion on a triangular lattice. 

We made an additional run of 10 8 samples and com- 
pared the density along the horizontal centerlinc (y = 
1/2), and found agreement with J7J within 1%, with 
slightly larger deviations at the two endpoints (bound- 
aries). We find similar agreement of the simulated one 
point density with ®. 

Thus the factorization in Q agrees with simulation 
results, supporting our arguments for discarding the so- 
lution of lO giving rise to the 0x,5 term. We also find 
numerically that this factorization holds asymptotically 
(with the same value of C) when one or both anchor 
points is not on the edge. (In this case the change of C 
from 1 to the quoted value occurs over a number of lattice 
spacings that scales with lattice size, indicating that the 
rhs of (Q has an additional term.) Furthermore, CFT 
calculations indicate that a similar factorization applies 
to the critical Potts models as well. Hence these results 
seem to apply more generally. 

Note that CFT implies that the factorization in (J7J 
remains valid under any conformal mapping. 

To understand these results better we combine and 
(|7|). which, in an obvious notation, gives 

P2 P ts{z;x a ,x b ) = j^^p p t{z;x a )p pt (z-,x b ). (8) 

Now the densities in JSJ, (and the factor 1/D 1 / 3 ), may be 
understood as the probability V of a cluster that connects 




FIG. 1: (color online) Simulation results for p p t(z;y) and 
P2 P ts(z;y a ,yb). 



the points in question (or the anchor points). Now since, 
as mentioned, "0 reduces to on a free boundary, if we 
take z — > x c , JEl reduces to 



(<t>(x a )4>(x b )(f)(x c )) = 

C\J (4>(x a )4>(xb)} (4>(x a )4>(x c )) (4>(xb)4>(x c )} , 



(9) 



where C is exactly the (boundary) OPE coefficient of the 
term that we retained in the solution of (J5J). Therefore 



V(z, x ai x b ) = C \JV{x ai x h )V{z,x a )V{z,x h ), (10) 

with C as above. Except for the square root, i|10|) resem- 
bles the Kirkwood superposition approximation familiar 
from the theory of fluids [l8j, which has been applied to 
percolation However it should be emphasized that 
(II 01) is both exact and universal. 

Finally we consider a cluster anchored along the entire 
interval xi < x < X2- Arguing as above, the cluster den- 
sity is now given by the four-point correlation function 

Pint(z;X 1 ,X 2 ) = (4>l,2{xi)(f>l,2(x2)lp(z)lp(z)) . (11) 

Proceeding as for |0J , and redefining 77 by replacing x a — > 
x\ and x b — > £2, gives 



(12) 
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F(rf) = 0. 



Since 12 is a level-two operator <|12fl is second order. 

The limit x%, X\ — > gives a single cluster anchored at 
the origin. The solution of (|12f> goes either as F_ ~ 77 



1/3 
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or as F + ~ 77° in this limit, corresponding, respectively, 
to the appearance of <fr or of the unit operator in the 
OPE, as mentioned above. Retaining _F_ only, we find 
that (fHE|) (see below) reproduces QJ, as it should. 

Note that setting = reduces <|12|) to Cardy's dif- 
ferential equation for the crossing probability |lj, as it 
should, since that quantity involves four boundary oper- 
ators 0i,2 with dimension zero. 

The two independent solutions F + and F_ of (|12|l may 
be written as 



*±(»7) = 
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We next set z = x + iy and X\ = —L/2,x 2 = L/2 so 
that the cluster is anchored on an interval of length L 
centered at the origin. Substituting these values into the 
cross-ratio 77, one finds, with appropriate choice of the 
square root, 



Pint(z; L) 



1 



,5/48 



(14) 



L 2 



4(x 2 + y 2 ) 



v/L 4 - 8L 2 (a; 2 - y 2 ) + I6(x 2 + y 2 ) 2 J 

Note that for L -> 0, p lnt (z;L) -> (L 1 / 3 /2 1 / 6 ) p pt (z), 
as expected since the interval reduces to a point. 

In the limit y — ► 0, (|14f> becomes a power series in 
y 2 when |i| < L/2 and y 1 / 3 times a power series in y 2 
when I a; I > L/2. This shows that in the boundary OPE 
[23,1^1 for tp, to leading order an operator of dimension 
h = (the unit operator, presumably) appears for fixed 
boundary conditions while an operator of dimension h = 



1/3 {4>, by the result above) appears for free boundary 
conditions. The latter conclusion emphasizes again that 
in order for ip to measure the density it must connect 
to a cluster. Note that an analogous dependence of the 
boundary OPE on boundary conditions occurs for the 
magnetization operator (</>2 2 or <f>i 2) in the Ising model 

m 

In order to make contact with 2-D electrostatics again, 
we first integrate $dip over the position of the dipole. Ex- 

plicitly, let ty(z;L) = J_ L / 2 ^dip(z; xo)dxQ. Note that if 
<&dip is taken as the imaginary potential as above, then 
^{z; L) is the imaginary potential for a finite dipole from 
—L/2 to L/2 (in the other interpretation it is the poten- 
tial of a dipole layer between these points). It is straight- 
forward to show that (|14fl may be expressed as 



Pint(z; L) 



y 



5/48 



,1/3 



(15) 



Again we see in (|15fl a kind of exponential additiveness, 
via the sine function. 

In summary, we present new results for the density 
of anchored clusters in two-dimensional systems at the 
percolation point. Our findings deepen and extend the 
known connections between the percolation point and 2- 
D electrostatics, and appear to be generalizablc. 
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